The role of differential rotation in the evolution of the r-mode instability 
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We discuss the role of differential rotation in the evolution of the I = 2 r-mode instability of 
a newly born, hot, rapidly-rotating neutron star. It is shown that the amplitude of the r-mode 
saturates in a natural way at a value that depends on the amount of differential rotation at the time 
the instability becomes active. It is also shown that, independently of the saturation amplitude of 
the mode, the star spins down to a rotation rate that is comparable to the inferred initial rotation 
rates of the fastest pulsars associated with supernova remnants. 
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I. Introduction 

The r-modes are non-radial pulsations modes of rotat- 
ing stars that have the Coriolis force as their restoring 
force and a characteristic frequency comparable to the 
rotation speed of the star 0. These modes are driven 
unstable by gravitational radiation in all rotating stars 
. A deeper understanding of r-modes and their astro- 
physical relevance requires taking into account nonlinear 
effects in the evolution of the r-mode instability. One 
such a nonlinear effect, differential rotation induced by 
r-modes, has been studied by several authors. Rezzolla, 
Lamb and Shapiro Q were the first to suggest that dif- 
ferential rotation drifts of kinematical nature could be 
induced by r-modes, deriving an approximate analytical 
expression for this differential rotation. Afterwards, the 
existence of such drifts was confirmed numerically both in 
general relativistic Q and Newtonian hydrodynamics . 
Recently, an exact r-mode solution, describing differen- 
tial rotation of pure kinematic nature that produces large 
scale drifts along stellar latitudes, was found within the 
nonlinear theory up to second order in the mode's ampli- 
tude In this paper we discuss the role of differential 
rotation in the evolution of the I — 2 r-mode instability 
of a newly born, hot, rapidly-rotating neutron star. 



II. The evolution model of Owen et al. 

A few years ago, Owen et al. Q proposed a simple 
model to study the evolution of the r-mode instability in 
newly born, hot, rapidly-rotating neutron stars. Within 
this model, it is assumed that the time evolution of the 
system (star and r-mode perturbation) is characterized 
by two parameters: the angular velocity of the star, Q{t), 
and the amplitude of the mode, a{t). The total angular 
momentum of the star is then given by 



J = m + Jc 



(1) 



where the momentum of inertia / of the equilibrium con- 
figuration is / = IMR^, with / = 0.261, and the canoni- 
cal angular momentum Jc of the r-mode perturbation is 
Jc = ~3/2a^njMR'^, with J = 1.635 x 10"^. Here, only 
the / = 2 r-mode is being considered and it is assumed 
that the mass density p and the pressure p of the fluid are 
related by a polytropic equation of state p = kp^, with k 
such that the mass and radius of the star take the values 
M = 1.4Mq and R = 12.53 km, respectively. 

Within this model it is also assumed that the total an- 
gular momentum of the star decreases due to the emission 
of gravitational radiation and that the physical energy of 
the r-mode perturbation (in the co-rotating frame) in- 
creases due to the emission of gravitational radiation and 
decreases due to the dissipative effect of viscosity. These 
assumptions lead then to a system of differential equa- 
tions for the angular momentum of the star, f2, and the 
amplitude of the r-mode, a: 



da 



2n Qo? 
n7l + Qa2' 

a. a 1 — Qc? 
TGR TV 1 + Qd^ ' 



(2) 
(3) 



where Q — 3J/(2/) = 0.094, the gravitational-radiation 
and viscous timescales, and ry, are given by 
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In the above expressions, the fiducial timescales are 
fgfl, = -3.26 s, fs = 2.52 x 10*^ s and = 6.99 x lO^s. 

For a newly born, hot, rapidly-rotating neutron star 
there is an interval of relevant temperatures and angular 
velocities of the star for which the gravitational timescale 
is much smaller than the viscous timescale, tgr ^ Ty. 
Therefore, for this interval of temperatures and angu- 
lar velocities, we can neglect in the right-hand side of 
Eqs. l|2Jl and © the terms proportional to Ty^ and ob- 
tain the solution = fio and a ~ exp{— (t — ^o)/''■Gi^}• 
If the initial angular velocity is chosen to be Qq = 
where fix — {2/3)^/TrGp is the Keplerian angular veloc- 
ity at which the star starts shedding mass at the equa- 
tor, then Tgr = —37.1s, implying that the perturbation 
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grows exponentially from an initial amplitude On = 10~^ 
to values of order unity in just about 500 s After 
this first stage of evolution of the r-mode instability, the 
amplitude a has to be forced, by hand, to take a cer- 
tain saturation value a^at ^ = 3.26; in the second 
stage of the evolution it is then assumed that a — asat 
and Q is determined from the equation 

dfl _ 2n 

dt TGR 1 

which yields the solution 



(6) 
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where is the time at which occurs the transition from 
the first to the second stage of evolution. During the 
second stage of evolution, the star loses its angular mo- 
mentum, spinning down to its final angular velocity. As 
can be seen from the above equations, the final angular 
velocity of the star depends critically on the saturation 
value of the mode's amplitude asat] for instance, after one 
year of evolution, for Sip = and asat = 1 one obtains 
n ~ O.lil/f, while for a^at = 10~^ the angular velocity 
is ~ Q.^VIk- The fact that the growth of the mode's 
amplitude has to be stopped by hand at a certain satu- 
ration value introduces an element of arbitrariness into 
the solution, permitting, for instance, that agreement be- 
tween the predicted final value of the angular velocity of 
the star and the value inferred from astronomical obser- 
vations can always be achieved by fine-tuning the value 
of the saturation amplitude. As will be seen in the next 
section, this arbitrariness disappears when differential ro- 
tation is taken into account. 



In the previous expression, K is & constant fixed by the 
choice of initial data and gives the initial amount of differ- 
ential rotation associated with the r-mode. The specific 
case K = —2, for which the physical angular momentum 
of the r-mode perturbation coincides with the canonical 
angular momentum, corresponds to the model of Owen et 
al. (]j described in great detail in the previous section. 
There is not, to our knowledge, any physical condition 
that forces K to take such a particular value K — —2. 
Therefore, we will consider in this paper arbitrary values 
of K in the interval -5/4 ^ if < 10" 

Following the procedure described in the previous sec- 
tion, we arrive at a system of two first-order coupled dif- 
ferential equations determining the time evolution of the 
amplitude of the r-mode a{t) and of the angular velocity 
of the star Q.{t): 



'dt 
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(10) 
(11) 



valid when the damping effect of viscosity can be ne- 
glected relatively to the driving effect of gravitational 
radiation. 

This system of equations was solved analytically and 
discussed in great detail in Ref. [9]. In the initial stages 
of the evolution of the r-mode instability a increases ex- 
ponentially 0, 



a{t) ~ aoexp<( 0.027 ( ^ 



sec 



(12) 



while for later times, a increases very slowly as 



III. Evolution model with differential rotation 

Differential rotation induced by r-modes does con- 
tribute to the physical angular momentum of the r-mode 
perturbation [g. Therefore, a model for the evolution of 
the r-mode instability should include the effect of differ- 
ential rotation. Here, we adopt the model of Owen et al. 
, described in the previous section, with the important 
difference^ that the total angular momentum of the star 
is given by 



(8) 



where the physical angular momentum of the r-mode per- 
turbation is 



'-a^n{4:K + b)jMR^ 



(9) 



a{t) ~ 2.48 



fix J V sec 



(13) 



The amplitude of the mode saturates in a natural way 
and a smooth transition between the regimes H12|) and 
(|13|l occurs for t—to ~ fewx 10^ seconds (see Fig.^. This 
contrasts with the model of Owen et al. , in which, af- 
ter the short initial period of exponential growth, the 
amplitude a has to be forced by hand to take a certain 
saturation value asat ^ Q~^^^. As can be seen from 
Eq. (fT!^ . the saturation value of the amplitude of the 
mode depends crucially on the parameter K, namely, 
asat (X {K + 2)^^/^. If if ~ 0, corresponding to a situa- 
tion in which the initial amount of differential rotation is 



^ There is another difference, albeit less significant, between the 
model of Owen et al. and the model we adopt here, namely, 
we deduce the evolution equations just from angular momentum 
considerations. 



^ The upper limit for K results from the fact that one wishes to im- 
pose the condition that \S^^^J{to)\ <C I^o (the initial amplitude 
CfQ = a(to) is considered throughout the paper to be 10~^). The 
lower limit for K results from the fact that we do not want to 
saturate the amplitude of the mode by hand, a procedure needed 
for the case K < —5/4 in order to avoid that the total angular 
momentum of the star becomes negative. 
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FIG. 1: Time evolution of the amplitude of the r-mode for 
different values of K. The initial values of the amplitude of the 
mode and of the angular velocity of the star are, respectively, 
Qfo — 10^^ and flo — fix- 
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FIG. 2: Time evolution of the angular velocity of the star for 
different values of K. The initial values of the amplitude of the 
mode and of the angular velocity of the star are, respectively, 
Qfo ~ 10~^ and flo = ^^k- 



small, then the r-mode saturates at values of order unity. 
If, on the other hand, the initial differential rotation as- 
sociated to r-modes is significant, then the saturation 
amplitude asat can be as small as 10"'^ — 10~^. 

Let us now turn our attention to the time evolution of 
the angular velocity of the star, H.. In the initial stages 
of the evolution of the r-mode instability il decreases as 



while for later times decreases slowly as 




The smooth transition between the regimes H14|l and 
((TK|l occurs for t — to ~ few x 10^ seconds (see Fig. |21l. 
Remarkably, in the later phase of the evolution, the an- 
gular velocity does not depend on the value of K and, 
consequently, does not depend on the saturation value of 
a. This contrasts with the results obtained in Ref. [7], 
where the value of O depends critically on the choice of 

asat- 

After about one year of evolution, when the dissipative 
effect of viscosity becomes dominant and starts damp- 
ing the mode, the angular velocity of the star becomes 
^^one year — 0.042fiif (for flo = flR), in good agreement 
with the inferred initial angular velocity of the fastest 
pulsars associated with supernova remnants. 



ity of a newly born, hot, rapidly-rotating neutron star. 
We have shown that, a few hundred seconds after the 
mode instability sets in, the amplitude of the r-mode 
saturates in a natural way at values that depend on the 
initial amount of differential rotation associated to the 
r-mode. If the initial differential rotation of r-modes 
is small, then the r-mode saturates at values of order 
unity. On the other hand, if the initial differential rota- 
tion is significant, then the saturation amplitude can be 
as small as 10~^ — 10~^. These low values for the sat- 
uration amplitude of r-modes are of the same order of 
magnitude as the ones obtained in recent investigations 
on wind-up of ma gne tic fields j^f and on nonlinear mode- 
mode interaction • We have also shown that the value 
of the angular velocity of the star becomes, after a short 
period of evolution of the r-mode instability, very insen- 
sitive to the saturation value of the mode's amplitude. 
After about one year of evolution the angular velocity 
is 0. 042^1 (for any asat), in good agreement with the 
inferred initial angular velocity of the fastest pulsars as- 
sociated with supernova remnants. 
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